Single-spin magnet ometry with multi-pulse sensing sequences 
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We experimentally demonstrate single-spin magnetometry with multi-pulse sensing sequences. 
The use of multi-pulse sequences can greatly increase the sensing time per measurement shot, re- 
sulting in enhanced ac magnetic field sensitivity. We theoretically derive and experimentally verify 
the optimal number of sensing cycles, for which the effects of decoherence and increased sensing time 
are balanced. We perform these experiments for oscillating magnetic fields with fixed phase as well 
as for fields with random phase. Finally, by varying the phase and frequency of the ac magnetic field, 
we measure the full frequency-filtering characteristics of different multi-pulse schemes and discuss 
their use in magnetometry applications. 
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The ability to sense weak magnetic fields with nanome- 
ter scale resolution has important applications in funda- 
mental and biomedical sciences as well as information 
storage technology. Several architectures for highly sen- 
sitive magnetometers have been implemented, such as su- 
perconducting quantum interference devices (SQUIDS) 
|H , Hall sensors [1| , sensors based on magnetic resonance 
force microscopy (MRFM) [3[ and atomic vapors [4] . Re- 
cently, approaches to magnetometry based on tracking 
the evolution of a single electron spin have been proposed 
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The magnetic field to be detected shifts the energy 
levels of the spin through the Zeeman effect. When a 
superposition of spin states is prepared, the Zeeman shift 
leads to a phase difference proportional to the magnetic 
field, which can be detected as a population difference 
after application of a suitable control pulse. 

The nitrogen vacancy (NV) defect center in diamond 
is a prime candidate for single-spin magnetometry [6|- 
0]. The NV center combines an excellent field sensitiv- 
ity with high spatial resolution resulting from the near- 
atomic size of the defect. It has a ground state para- 
magnetic spin (5=1) with the m s = and m s = ±1 
levels split by 2.87 GHz at zero magnetic field, with 
the quantization axis oriented along the symmetry axis 
of the defect. The NV center spin can be initialized 
and its spin population difference can be detected op- 
tically by measuring its spin-dependent photolumines- 
cence (PL) [To[. First experimental demonstrations of 
NV-based magnetometry using a sensing sequence with 
a single 7r-pulse (spin echo) have outlined its potential as 
an ultra-sensitive detector [8|; even better performance 
can be achieved by applying more control pulses @, [H| . 

Here, we present a detailed experimental study of sin- 
gle NV center magnetometry with multi-pulse sensing 
sequences. We demonstrate that multi-pulse sequences 
greatly improve the sensitivity to oscillating magnetic 
fields. We find the optimal number of control pulses 



( a ) control 




(b) 

t* ° 

+1 



BM) 



4fc - 1 T 



+ Y(t) 



B z (t) 01 



toggling 
►frame 



FIG. 1: (a) Magnetometry setup. A CPW transmission line, 
fabricated directly on a lb diamond sample, is connected to 
control lines on both ends. One control line supplies 2.5 GHz 
microwave bursts for spin control. For details of this part of 
the setup see Ref. [l3[. We use the other side of the CPW 
to supply the low frequency (< 5 MHz) ac field B z (t) that 
we aim to detect. A single NV center located in between the 
conductors of the CPW is used as the magnetometer, (b) 
Multipulse magnetometry. An ac field B z (t) modulates the 
phase of the NV electron spin. Timing 7r-pulses such that 
they coincide with the nodes of B z (t) effectively multiplies 
the field by Y(t). The bottom trace shows the resulting B z (t) 
in the toggling reference frame, which flips whenever a 7r-pulse 
is applied to the spin. 



as function of field frequency and spin-echo decay time. 
These studies are performed both for fields with known 
phase and for fields with a phase that randomly fluctu- 
ates between measurement shots. We finally show that 
multi-pulse sequences can be used to achieve a high de- 
gree of frequency tunability and selectivity by exploiting 
their frequency-domain characteristics. 

We use a single NV center in a nitrogen-rich lb bulk di- 
amond sample (Element Six) as our magnetometer. The 
setup is shown schematically in Fig. 1(a). We apply 
a static magnetic field oriented along the NV center's 
quantization axis z that allows us to selectively address 
the m s = <rt m s = — 1 spin transition. Within this 
subspace, the NV spin is equivalent to a spin- 1/2. Us- 
ing lithographically defined on-chip coplanar waveguides 
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(CPWs) [12] we achieve high fidelity («99%) control of 
the NV center electron spin state 13|, |l4[ . The length of 
a single 7r-pulse used here is t p = 8 ns. With the CPW 
low-frequency magnetic fields can be applied along the 
z-axis as well. 

We consider a time- varying oscillating magnetic field 
oriented along the quantization axis z of the probing spin 
B z (t) — b z sm(27rft + </>). Components oscillating along 
the transverse axes will average to zero due to the rapid 
(2.5 GHz) precession of the NV spin. The working prin- 
ciple behind detecting a field of the form B z (t) is outlined 
in Fig. 1(b). An electron spin initialized in a superpo- 
sition state precesses under the influence of the oscillat- 
ing field. During each half-cycle the electron spin phase 
acquires 5§ = J^ 2f 2ir^B z (t)dt = (27^//) cos 0, with 
7 = 28 GHz/T for the electron spin of an NV center. 
Since the field oscillates the total phase after many cy- 
cles averages to zero. This can be prevented by using an 
TV-pulse sequence with evenly spaced 7r-pulses 



[t/2 - 7T - T - 7T - T/2] 
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Note that this sequence has the same timing properties as 
the Carr-Purcell sequence (CP) [IH known from NMR. 
CP-like sequences have very recently been explored with 
NV centers in the context of dynamical decoupling from 
decoherence by a spin bath environment [U], Q, fl6| . 
When the 7r-pulses of the sequence coincide with the 
nodes of B z (t), the electron accumulates a phase which 
increases with the length of the sequence. This can be 
understood by moving to the toggling reference frame of 
the electron spin from where B z (t) — >> —B z (t) after each 
7r-pulse. Formally, B z (t) is multiplied by a time-domain 
filter function Y(t) which changes sign each time a tt- 
pulse is applied (Fig. 1(b)). After N 7r-pulses (or half- 
cycles) the total phase of the electron spin state becomes 
A$ = N6&. 

We first analyze the case when B z (t) always has the 
same phase relation with respect to the sequence of ir- 
pulses (0 = 0) @, 0. This corresponds to a situation 
where the phase of the field to be measured is under con- 
trol of the experimenter (e.g. when one aims to detect 
spins that can be adiabatically inverted periodically). In 
this case, for a spin initialized along the x-axis, the read 
out can be performed by rotating the final state by tt/2 
around the same axis (see Fig. 2(a)). The resulting signal 
is then S z (b z ) = \ sin [(2Njb z / f) cos 4>] after normaliza- 
tion of the PL levels to [— |, |] . 

In Fig. 2(b) we monitor the evolution of the NV elec- 
tron spin under application of an ac field with frequency 
/ = 1 MHz for increasing TV. The ac field is phase- 
locked such that the 7r-pulses coincide with the nodes of 
B z (t). We observe oscillations which demonstrate spin 
precession in the applied ac field. The amplitude of the 
oscillations decays exponentially due to decoherence and 
pulse imperfections. We can limit the influence of pulse 
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FIG. 2: Phase- locked magnetometry. (a) Measurement 
scheme. A spin initialized along the x-axis will accumulate 
a phase during the sensing stage due to B z (t). A final rota- 
tion around the x-axis is applied to transform the phase into 
a population difference, (b) A 1 MHz ac field with constant 
amplitude is measured with increasing number of pulses N. 
Solid line is a fit to Eq. ([2j), yielding b z = 1.6 fiT. (c) Signal 
for three ac fields with different frequency. The amplitude 
is rescaled by 1/N. Solid lines are fits to Eq. ([2]) yielding 
2Njb z /f « 0.2tt and T 2 = (2.86 ± 0.04) /is (d) Sensitivities 
calculated from the data (points) and from the fits (solid lines) 
in (c). 



imperfections on the signal decay by using the XY4 se- 
quence which is self-compensating for pulse errors [Tt( | 
and has the same timings as Eq. ([1]). We find the maxi- 
mum number of pulses that can be applied before pulse 
errors start to play a role to be ~ 130 pulses 

In our type lb diamond sample the coherence time is 
limited by dipolar interactions with electron spins, which 
are located at the sites of substitutional nitrogen atoms 
in the diamond lattice [18]. The one-pulse spin-echo sig- 
nal decays as - exp[-(r/T 2 ) 3 ] where T 2 = (2.8 ± 0.1) 
[is for the NV center used here [14]. Increasing the 
number of pulses to N reduces the signal by a factor 
~ exp[— 7V(t/T2) 3 ] for the CP sequence, which leads to 
the observed exponential decay in the signal in Fig. 2(b). 
For B z (t) oscillating with frequency / = [2(r + t p )] _1 ~ 
(2r) _1 (assuming t p <C r) this gives a total signal of 
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The optical detection of the spin population is limited by 
shot noise, which depends on the experimental parame- 
ters such as the number of photons collected per measure- 
ment shot $ and the contrast C which combine to give 
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the noise per measurement shot as z ~ l/(Cy/^) fA fl9|. 
The sensitivity for detecting a field oscillating in phase 
((f) = 0) is given by combining the shot-noise limited min- 
imum detectable field b m i n = as z \db z /dS z \ « as z ^ an d 
the total integration time T = N(r + t p ) = N/2f to give 



7l(f,N)=b min VT: 
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It is instructive to consider two limiting cases. If 
N <C 2/T2, decoherence is negligible and the use of a 
multi-pulse sequence improves the sensitivity by a factor 
1/y/N. In the other extreme, where TV ^> 2/T2, de- 
coherence has a detrimental influence on the sensitivity: 
r] oc e N /( 2 f T2 ^ . Thus, for a given frequency of B z (t) there 
exist an optimum number of pulses iVo P t . By minimizing 

77wefindJV pt=4T 2 3 / 3 . 

In Fig. 2(c) we demonstrate the detection of ac mag- 
netic fields of three different frequencies, for increasing 
N. In order to keep the signal in the linear regime the 
amplitude of the field is rescaled as N increases so that 
the total acquired phase after applying the sequence re- 
mains constant. For every measurement shot we also 
measure the zero-field signal to account for possible drifts 
in the setup. From the curves in Fig. 2(c) the sensitiv- 
ity is calculated and depicted in Fig. 2(d), along with the 
predicted iVo P t (calculated using the independently deter- 
mined T2 = (2.8 ± 0.1)/is). We observe that there indeed 
exists an optimum number of pulses for each frequency 
and find excellent agreement between the predicted N opt 
and the data. 

We now turn to the case where B z (t) has a phase 
which varies randomly between measurement shots, but 
remains constant during each individual measurement 
shot. The phase that a spin, initialized along x, acquires 
will be different for each measurement shot with zero av- 
erage, as depicted in Fig. 3(a). Rotating the final state 
over the x axis therefore yields zero signal when aver- 
aged over many measurements. However, as shown in 
Fig. 3(a), the state can also be rotated around the (or- 
thogonal) y-axis. This will give a signal, when averaged 
over the phase, of 
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where ()^ denotes averaging over cj> and Jo is the zeroth- 
order Bessel function of the first kind. In the absence of 
any field S z (0) = \. The data in Fig. 3(b) demonstrate 
that we can also detect the ac magnetic field if it has a 
random phase. 

A consequence of measuring the x-projection is that 
\dS z /db z \ « 2N 2 j 2 b z /f 2 will vanish as b z -> 0, leading 
to a divergence in the differential sensitivity. We there- 
fore turn to calculating the interval sensitivity [20| by 
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FIG. 3: Magnetometry for ac fields with random phase, (a) 
Measurement scheme. In each measurement shot the spin, ini- 
tialized along the x-axis, accumulates a different phase during 
the sensing stage that depends on f A signal can nonethe- 
less be measured by applying the final n/2 rotation over the 
y-axis. (b) Measured signal for a field with constant am- 
plitude and averaged over cj) for increasing N. Solid line is 
a fit to Eq. (|4}; the expected signal for — (dashed) is 
included as a reference, (c) Signal intervals for a random- 
phase ac field of 651 kHz. The amplitude b z is rescaled by 
1/N. Solid line is a fit to Eq. (0} yielding 2N~fb z /f « 0.47tt 
and T2 = (2.77±0.05)/is. (d) Interval sensitivities calculated 
from the data (points) and from the fit (solid line) in (c). 



extending the definition of the minimum detectable field 
per measurement shot to that measured for a given sig- 
nal interval: b m i n = (Js z b z / \S z (b z ) — 5^(0)|. The interval 
sensitivity rji for measuring oscillating fields with random 
phase is given by 



I Jo (2Njb z /f) 




(5) 



Analogously to the phase- locked experiment, we verify 
this expression experimentally in Fig. 3(c). The corre- 
sponding calculated interval sensitivities are depicted in 
Fig. 3(d). It shows a qualitatively similar picture as for 
the phase- locked case. Since we rescale b z by 1/N, also 
here N opt = 4T| f 3 . Again, we observe that the theory 
gives an excellent description of the data. 

Until now we discussed the situation where the se- 
quence was tuned exactly in resonance with B z (t). In 
order to analyze what happens when the sequence is de- 
tuned with respect to B z (t) (so r / ^) we move to 
the frequency domain. For a sequence of evenly spaced 
pulses and = the response of the magnetometer in 
the frequency domain is given by Q 
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FIG. 4: Frequency domain analysis, (a) Measured filter re- 
sponse Yn(/, t) of the CP sequence vs. frequency for a field 
amplitude of b z « 0.8 fiT for = and r = 192 ns. The 
scheme for measuring a phase- locked signal from Fig. 2(a) is 
used. A maximum signal is observed at / = 2(r + t p ) _1 = 
2.5MHz. As N is increased the bandwidth reduces and peak 
signal becomes higher. (b) Line traces from (a). Solid 
lines are fits using Eq. J6}- (c) Magnetometer response vs. 
frequency and phase using the CP-like XY4 sequence with 
N = 20 pulses and 8.16 fis integration time (r = 400 ns). 
The field amplitude is adjusted to ~ 1.7 fiT. The scheme for 
measuring a random phase signal from Fig. 3(a) is used, (d) 
Signal for the same field and integration time measured in 
(c) but now using the 20-pulse aperiodic UDD sequence. The 
upper panels in (c) and (d) depict the signal averaged over 
the phase. 



Figure 4(a) depicts the signal detected as a function 
of frequency of B z (t) (with cj) = 0) for N = 4 to 80 for 
fixed r, mapping out the complete filter function Eq. (j6]). 
Figure 4(b) shows line traces for three different N . With 
increasing N, the bandwidth of the response decreases 
by a factor N while the peak signal increases N times. 
Therefore, by tuning r a single frequency component can 
be selected. This is useful for measuring the linewidth of 
B z (t) or for spectroscopic applications. The resolution is 
set by the full-width at half-maximum (FWHM) A/ « 
0.3/(7Vr). 

We study the influence of the phase of B z (t) in more 
detail by performing a measurement at a fixed number 
of pulses for a range of initial phases. The results are 
depicted in Fig. 4(c). By averaging over the phase the 
response to fields with random phase is retrieved (see top 
panel of Fig. 4(c)). The filter function of sequences other 
than CP-like sequences can be investigated in a similar 
way. 

The Uhrig dynamical decoupling (UDD) sequence 21] 



has been conjectured as a valuable tool in detecting ran- 
domly fluctuating fields [22]. The TV-pulse UDD se- 
quence has pulses spaced irregularly according to = 
sin 2 [7r/c/(27V + 2)] with the fc-th pulse spacing. A char- 
acterization of the UDD sequence similar to Fig. 4(c) is 
presented in Fig. 4(d). UDD is seen to give a broader 
frequency response and reduced peak signal compared to 
a sequence with CP-like timings (Eq. (pQ)). The higher 
peak signal and reduced bandwidth of the latter is es- 
pecially useful when gradients are used to achieve high 
spatial resolution. It will in general depend on the nature 
of B z (t) and the specific application for which sequence 
the best performance will be achieved. 

In conclusion, we have reported a detailed investigation 
of spin-based magnetometry with multi-pulse schemes. 
Our results show significantly enhanced performance 
both for ac fields with known and with unknown phase. 
These results pave the way towards unprecedented mag- 
netic field sensitivity beyond the limit set by the spin- 
echo sequence. Note that the multi-pulse sequences also 
make the magnetometer insensitive to instabilities in the 
setup, such as drifts in the applied static magnetic field 
or in temperature [23] . The insights gained here will help 
guide experimenters in tailoring the pulse sequence and 
number of control pulses to their specific application. 
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